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Abstract 

In superstring theories, there exist various dilaton and modulus fields which 
masses are expected to be of the order of the gravitino mass m 3 / 2 - These 
fields lead to serious cosmological difficulties, so called "cosmological moduli 
problem", because a large number of moduli particles are produced as the 
coherent oscillations after the primordial inflation. We make a comprehen- 
sive study whether the thermal inflation can solve the cosmological moduli 
problem in the whole modulus mass region ~ 10 eV-10 GeV predicted 
by both hidden sector supersymmetry (SUSY) breaking and gauge-mediated 
SUSY breaking models. In particular, we take into account the primordial in- 
flation model whose reheating temperature is so low that its reheating process 
finishes after the thermal inflation ends. We find that the above mass region 
m^,(~ 1^3/2) ~ 10 eV-10 GeV survives from various cosmological constraints 
in the presence of the thermal inflation. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is one of the most attractive extensions of the standard model. 
In virtue of the SUSY the electroweak scale can be stabilized against the radiative correc- 
tions. Furthermore, in the SUSY grand unified theories the unification of the standard gauge 
couplings can be realized. 

However, overviewing the cosmology of the SUSY model, we are faced with various 
difficulties. One of the cosmological problems is the gravitino problem PHI. This problem 
still exists even if the universe experienced a primordial inflation, since the gravitino is 
reproduced by scatterings with particles in the thermal bath at the reheating epoch. In the 
hidden sector SUSY breaking (HSSB) models f| the gravitino has a mass of the order of the 
electroweak scale: m.3/2 ~ 10 2 -10 3 GeV, and it decays soon after the big bang nucleosynthesis 
(BBN). Since high energy photons produced by the gravitino decay might destroy the light 
elements (D, 3 He, 4 He) synthesized by the BBN, the reheating temperature of the primordial 
inflation should be low enough not to conflict with the observations (e.g., see a recent 
analysis in Ref. ||). On the other hand, the gauge mediated SUSY breaking (GMSB) 
models predict a light gravitino of mass m 3 / 2 ~ 10 eV-1 GeV. If the gravitino mass is 
m 3 / 2 > IkeV jLf, a low enough reheating temperature of the inflation is also required to 
avoid the overclosure by the stable gravitino 0. 

Furthermore, when one considers the SUSY models in the framework of the superstring 
theories they suffer from a more serious problem, i.e., "cosmological moduli problem" 
P~p|. As a general consequence of the superstring theories, various dilaton and modulus 



fields appear. (We call them as "moduli" throughout this paper.) These moduli fields 
cf) are expected to acquire masses of the order of the gravitino mass (m^ ~ ^3/2) from 
nonperturbative effects of the SUSY breaking fLl| . Their lifetimes, since they have only the 
gravitationally suppressed interaction, are roughly estimated as 

~ N ~ lM t - !0 17 «* AT- M^— V 3 , (1) 
m\ V100 MeV/ ' v ' 

where N denotes the number of the decay channels and M v \ is the Planck scale M v \ = 
1.2 x 10 19 GeV. The existence of moduli <fi with such long lifetimes leads to various cosmo- 
logical difficulties, because a large number of moduli particles are produced as the coherent 
oscillations after the primordial inflation. In the HSSB models decay soon after the BBN 
as the gravitino, and hence the light elements might also be destroyed pHll"|j. On the other 



hand, for the lighter stable moduli predicted by the GMSB models, say < 100 MeV, the 
energy of the oscillation lasting until the present overdoses the universe. Moreover, moduli 
with masses m$ ~ 0.1 MeV-1 GeV give too much contributions to the x(7)-ray background 



spectrum JT2|. Therefore, moduli particles typically predicted by both HSSB and GMSB 
models bring a cosmological disaster. 

One way to evade these difficulties due to the string moduli is to give heavy masses 
(<p ~ 100 TeV) to all moduli, so that the decays of moduli take place before the BBN fl3 



and the lightest superparticles produced by the moduli decays do not overdose the universe 
|14| . Although a large entropy is produced by the moduli decay, the present observed 



baryon asymmetry can be naturally explained by the Affleck-Dine mechanism |T5[ as shown 
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in Ref. fl6f . However, such heavy moduli masses (i.e., a heavy gravitino mass) could be only 
achieved by some specific models. 

Here we would like to emphasize that this cosmological moduli problem could not be 
solved by the primordial inflation, even if one assumed an extremely low reheating temper- 
ature about 10 MeV which is limited from below by the BBN observation. This is a crucial 
difference from the gravitino problem. Therefore, we require some extra mechanism other 
than the primordial inflation to dilute the moduli mass density sufficiently^ 

One of such dilution mechanisms is the thermal inflation model proposed by Lyth and 
Stewart \TE\^ The thermal inflation occurs before the electroweak phase transition and 



produces tremendous entropy (just) before the BBN epoch, which leads to sufficient dilution 
of the moduli density. In Ref. |18| it was shown that the moduli problem in the HSSB models 
can be solved by the thermal inflation.^ 

On the other hand, the attempts to solve the problem in the GMSB models by the 
thermal inflation were done in Refs. |T^|T^J2^J2^] . It was shown by Ref. [TP] that it can 



solve the overclosure problem for the stable moduli whose masses are ~ 10 eV-10 1 



GeV. However, in Refs. [^,^,^], the stringent x(7)-ray background constraint was found 
to exclude the modulus mass region ~ 10 _4 -1 GeV. Thus the thermal inflation can 
cure a part of the cosmological difficulties of the string moduli predicted by both HSSB and 
GMSB models. 



However, most of previous works (except Ref. []T8j ) assumed that the reheating process 
of the primordial inflation be completed before the moduli oscillations start, i.e., only the 
inflation models with relatively high reheating temperatures were considered. Even in Ref. 



the reheating process is assumed to end before the thermal inflation starts. Without 
this restriction one might obtain wider allowed region for the modulus mass. In this paper, 
therefore, we make a comprehensive study whether the thermal inflation can solve cosmo- 
logical problem of the moduli particles with masses ~ 10 eV-10 4 GeV (i.e., the whole 
modulus mass region predicted by both HSSB and GMSB models), especially taking into 
account the inflation model whose reheating temperature is low enough so that its reheat- 
ing process finishes after the thermal inflation ends. We find that the above mass region 
m^(~ ^3/2) ~ 10 eV-10 4 GeV survives from various cosmological constraints^ if we con- 



1 Some other approach to the solution has been proposed IfTfll , but it is difficult to construct a 
realistic model. 



2 The other possible mechanism is the oscillating inflation by Moroi [20| which works in the GMSB 
models. A mini-inflation takes place while a scalar field corresponding to the flat direction oscillates 
along the logarithmic potential induced by the GMSB mechanism and it significantly dilutes the 
relic density of the moduli. See also Ref. 



3 More precisely, the thermal inflation can solve the moduli problem if one assumes the primordial 
inflation model which reheating process completes after the moduli begin the coherent oscillations. 

4 Here we do not include a constraint from the present baryon asymmetry which excludes some 
modulus mass regions. See the discussion below. 
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sider the primordial inflation with sufficiently low reheating temperature in addition to the 
thermal inflation. 

The organization of this paper is as follows. In Sec. [TT]we first review about the cosmolog- 
ical difficulties of the string moduli whose masses are ~ 10 eV-10 4 GeV. The (original) 
thermal inflation model proposed by Lyth and Stewart is explained in Sec. [TTI|. Then, in Sec. 
[V , we study whether this original thermal inflation could solve the moduli problem, con- 
sidering various models of the primordial inflation. In the original thermal inflation model, 
however, there appears an i?-axion which is the NG boson from the spontaneous breaking 
of the i?-symmetry. In the GMSB models the i?-axion is light enough so that the flaton, 
which causes the thermal inflation, almost decays into them. Then, as shown in Ref. |23fl . 
the original model could not produce a sufficient entropy to dilute the light moduli predicted 
by the GMSB. To avoid this difficulty we introduce in Sec. [V] the "modified" thermal in- 



flation by Ref. |23|, which forbids the flaton decay into i?-axions. Then, we investigate the 
solution of the moduli problem by this modified model in Sec. |VI[ Finally, in Sec. |V11| , we 
discuss two extra problems when one assumes the tremendous entropy production enough 
to dilute the string moduli, i.e., problems of baryogenesis and dark matter. We find that 
the problem of the baryon asymmetry is so serious that we have only two possible solution 



so far; (i) sufficient baryons can be generated through the Affleck-Dine mechanism for 
the modulus mass region m^(~ m 3 / 2 ) < 1 MeV and (ii) for m^,(~ ^3/2) ~ 1 10 4 GeV the 
electroweak baryogenesis can produce the present observed baryon asymmetry. 



II. COSMOLOGICAL DIFFICULTIES OF MODULI PARTICLES 

First of all, we explain the cosmological difficulties due to the string moduli particles. 
In the following analysis we assume only one modulus field, 0, with the mass ~ 7713/2 
to derive the conservative constraints, although we can easily extend the analysis for more 
general cases. 

After the primordial inflation ends, the modulus field <fi is considered to be displaced from 
the true minimum and the displacement is of the order of Mq [Mg is the reduced Planck 
scale M G = 2.4 x 10 18 GeV). This is because of the additional SUSY breaking effect due to 
large vacuum energy of the inflaton. Then, when the expansion rate of the universe (i.e., 
the Hubble parameter H) becomes comparable to the modulus mass, the modulus starts to 
oscillate around its true minimum with the initial amplitude <fto ~ Mq- At this time, the 
energy density of this oscillation is given by p^ = rrv^cf)^^ and the cosmic temperature of 
the universe is estimated as 

Tt K ( JL) ^ , ,2 x 1* GeV (^f , (2) 

where g*(— 200) counts the effective degrees of freedom of the radiation. Thus the ratio 
between p^ to the entropy density s is given by 

PJL s H o, 89 x 10 » GeV f J=U ( A V . (3) 

T?9.I? VI GeVT 1 \M C ) 

Here it should be noted that p^ and s are both diluted like R~ 3 (R is the scale factor of the 
universe) as the universe expands. Thus this ratio takes a constant value until today (or 
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until the modulus decays) if "no" extra entropy is produced. Since, at T = T^, the energy 
density of the radiation, p R , is comparable to that of the modulus (p R ~ T| ~ p^) and 
is diluted faster as R~ 4 , the modulus oscillation soon dominates the whole energy of the 
universe. 

In deriving Eq. (|3|) we have assumed that the reheating process of the primordial inflation 
be completed before the beginning of the modulus oscillation, i.e., the decay rate Y 9l of the 
inflaton be larger than m^.Q On the other hand, when T Vl < m<j>, the energy of the modulus 
oscillation is expected to be diluted by the entropy production of the primordial inflation. 
In this case, when the reheating process completes at H ~ r w , the ratio p^/s is estimated 
as 



P± 
s 




(4) 



Here Trj denotes the reheating temperature of the primordial inflation which is given by 

1/4 



RI 



90 



t vi m g . (5) 



In order to keep the success of the BBN it should be higher than about 10 MeV and we 
obtain 

& J> 1.25 x 10- GeV (A)'. (6) 

Comparing with Eq.(|3]), it is seen that the primordial inflation with a low Trj does dilute 
the modulus energy significantly. However, even if we assume the primordial inflation with 
extremely low reheating temperature, the energy of the modulus oscillation is still large. It 
leads to disastrous effect on the thermal history of the universe after the BBN if the modulus 
lifetime is long enough. 

The modulus, in fact, has a very long lifetime since it has only gravitationally suppressed 
interaction. From the naive dimensional analysis the modulus lifetime is estimated as shown 
in Eq. ([!]). Then we can roughly say that the modulus with a mass < 100 MeV has a 
longer lifetime than the present age of the universe ~ 10 17 sec. 

We can say more precisely about the modulus lifetime for the dilaton particle which is 
the most plausible candidate among various moduli. The dilaton has the following couplings 
to the kinetic terms of the gauge fields: 



AM, 



G 



5 In other words the reheating temperature of the primordial inflation is higher than given by 
Eq. (|). 

6 Note that the modulus oscillation is considered to start at least after the primordial inflation ends, 
since <C Hj where Hj denotes the Hubble parameter during the inflation and Hj = Vj / (\/3Mq) 
with the inflaton's vacuum energy V/. 
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where we introduce an order one parameter b which is determined by the string model 
and its compactification. For example, the dilaton has a coupling b = v2 |24| in some 
compactification of the M-theory ||25|| . If the dilaton mass is light enough, say < 1 GeV, 
it dominately decays into two photons ]l2j and the lifetime of the dilaton is estimated as 



r. = 7.6 x 10" 8 ec p . (8) 

Moreover, the dilaton could decay into two gluons similarly if kinematically allowed. Then 
the dilaton lifetime becomes shorter by a factor of 1/9 counting only the number of the 
final states. In addition, if is heavier than the electroweak scale, various decay modes 
may arise through the gravitational interaction with the SUSY standard model particles, 
and its decay lifetime becomes shorter. Comparing with the naive estimation Eq. ([!]), the 
both results are not different very much. Thus, in the present analysis we assume that the 
modulus lifetime be the same as the dilaton one and it decay into only two photons and two 
gluons if possible through the interaction (|7|) with 6 = 1 for simplicity. 

We turn to see constraints on the modulus energy from various cosmological observations. 
First of all, the stable modulus of mass < 100 MeV is constrained from the overclosure 
limit. The present energy density of the modulus oscillation should be smaller than the 
critical density of the universe p cr . This requirement leads to to 

El < ESL = 3.6 x lO" 9 h 2 GeV, (9) 
s s 

where s is the present entropy density and h denotes the present Hubble parameter in units 
of 100 km/sec-Mpc^ 1 . Note that even if the modulus is unstable, we obtain similar upper 
bound on the modulus abundance requiring that the decay products (i.e., photons or gluons) 
should not overdose the universe. Since the decay products are relativistic and its energy 
density decreases faster than the modulus one, the constraint on the abundance becomes 
weaker than Eq. (H). 

If the lifetime of the modulus is longer than the time of the recombination 10 12 
sec, i.e., < 1 GeV, a stringent constraint comes from the the observed cosmic x(7)- 
ray background spectrum |T2[]. The photon flux produced by the modulus decay directly 



contributes to the x(7)-ray backgrounds and we obtain the upper bound on the modulus 
abundance by requiring that they should not exceed the present observed spectrum. (One 
may find the detail analysis in Refs. |T2| , |23|| .) In fact this gives us a more stringent constraint 
than the above overclosure limit (§) for 1 GeV^ > 100 keV. Therefore, this is very 
stringent constraint on the models of the GMSB where ~ m%/2 ~ 1 GeV as first pointed 



out by Ref. [12 



The modulus abundance is also constrained from the spectrum of the cosmic microwave 
background radiation (CMBR). If the modulus lifetime is about 10 6 sec < < 10 12 sec, 
extra radiation energy produced by the modulus decay may cause the spectral distortion of 
the cosmic microwave background from the blackbody one. The observation by the COBE 



satellite p6| gives the following upper bounds on the modulus abundances: 

^ < 2.8 x lO- 3 GeV (^^y^ (10) 
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for 10 6 — - — sec < ta, < 1.4 x 10 9 — - — sec. 
V0.0125/ V - 0125 / 

^ < 2.6 x 10- 13 GeV (H^) 1/2 (11) 

for 1.4 x 10 9 ( j^g) sec < < 10 12 sec. 

Here Qb is the density parameter for baryons. 

Furthermore, the BBN plays a significant role to set limits to the moduli abundance. 
First, if the modulus oscillation exits at the cosmic temperature about 1 MeV, its extra 
energy accelerates the expansion of the universe and the weak interaction freezes out earlier. 
This leads to higher neutron-proton number ratio and overproduction of 4 He. To keep the 
success of the BBN, the modulus energy density should be smaller than that of one neutrino 
species [e.g., see Ref. M], i.e., 

^ < 1.2 x 1(T 4 GeV, (12) 
s 

for T/j, ^ 1 sec. 

Furthermore, when the modulus mass becomes larger than about 10 GeV, the modulus 
decays soon after the BBN. If the lifetime of the modulus is longer than about 10 4 sec, 
the high energy photons emitted from the modulus may destroy or overproduce the light 
elements of the universe synthesized by the BBN. Not to spoil the success of the BBN, the 
modulus abundance is stringently constrained [e.g., see the resent work ||]. Moreover, when 
the modulus lifetime becomes shorter than about 10 4 sec, the hadronic cascade processes 
associated with the modulus decay may modify the primordial abundances of light elements. 



This also puts the upper bound [27 



Although we have neglected the modulus decay into the SUSY standard model parti- 
cles, the stable lightest SUSY particle (LSP) produced by the modulus decay is potentially 
dangerous, since its energy density might overdose the universe. However, from the present 
mass limit of the LSP, we find that this constraint, if exists, is weaker than those from the 
distortion of the CMBR and the photo (hadronic) dissociation of the light elements. 

It should be noted that the modulus with ^ 10 TeV is cosmologically viable [|13 . 
From Eq. ([!]) the modulus with such a heavy mass decays before the BBN and reheats the 
universe as T > 10 MeV which gives the initial condition of the BBN. Then one can evade 
the previous constraints.)] 

We show these upper bounds on the modulus abundance Q<pf\ in Figs. |l], || and || taking 
the modulus coupling as b = 1, 0.1 and 10, respectively. The predicted modulus abundances 



7 However, in order to avoid the overclosure by the LSP produced by the modulus decay > 100 
TeV is required fl4|| . Furthermore, if the modulus decays into lighter gravitinos of mass 7713/2 ~ 100 



GeV-10 TeV, we are faced with the gravitino problem [28| 



8 For the stable modulus of mass < 100 MeV, $7^ is given by = (p^)o/p C r with the present 
energy density of the modulus (p</>)o- On the other hand, for the modulus of mass > 100 
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b=l 




[GeV] 

FIG. 1. Various cosmological upper bounds on the modulus abundance for the case that the 
modulus coupling is b = 1. The dotted line represents the upper bound from the overclosure limit 
on the abundance of the modulus (or its decay products) for < 200 GeV, and the upper bound 
from the BBN speed up effects for > 200 GeV. The dot-dashed line represents the upper bound 
from the x(7)-ray backgrounds. The short dashed line represents the upper bound from the CMBR 
spectrum. The long dashed line represents the upper bound from the dissociation of the BBN light 
elements. Note that there exist two typical modulus masses; ~ 100 MeV (the modulus lifetime 
is equal to the age of the universe) and m& — 1 GeV (the modulus decay into two gluons starts to 
open.) We also show the predicted modulus abundances of <j>Q = Mq by the solid lines for the case 
Triff, < T VI [Eq. (|3|)] and for the case > T VI with the reheating temperature Trj = 10 MeV [Eq. 
(§]• 

Eqs. (|3]) and (f|) with Trj = 10 MeV are also found in those figures. One can easily see that 
the string modulus with a mass from 10 eV to 10 TeV is excluded by the various cosmological 
observations. Since ~ m 3 / 2 , the whole gravitino mass region typically predicted by both 
GMSB and HSSB models is not cosmologically allowed. This difficulty is often referred as 
"cosmological moduli problem". Here we would like to stress that this problem could not 
be solved by choosing the model of the primordial inflation, i.e., even if one assumes the 
extremely low reheating temperature Trj ~ 10 MeV. This is very different from the gravitino 
problem. Therefore, we required some extra mechanism to dilute the modulus mass density 
sufficiently other than the primordial inflation. In the following, we consider the thermal 
inflation model proposed by Lyth and Stewart as such a dilution mechanism. 



MeV, fi^ is regarded as the ratio, (p<f> / s) d / (p C r / so) , where (p ( j ) /s)D denotes the ratio of the energy 
density of the modulus to the entropy density when the modulus decays. 
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FIG. 2. Same figure as Fig. [I] except for the value of the modulus coupling b = 0.1. 



b=10 



io z 



io 1 



io 1 



io 5 

a 

10° 



10 



10" 



-5 r 



fi.h(m 4 <r ) 



fi.h(m>r ) 



overdose 



X" 



/ 

/ /-I 



x(y)-ray 



CMBR.,^ 
r" - y BBN 



io" 8 io" 7 io" 6 io" 5 io" 4 io" 3 io" 2 10"' io° io 1 io 2 io 3 io 4 



[GeV] 

FIG. 3. Same figure as Fig.|] except for the value of the modulus coupling b = 10. 



9 



III. ORIGINAL THERMAL INFLATION MODEL 



In this section we review an original thermal inflation model proposed by Lyth and 
Stewart The thermal inflation is caused by a scalar field (called "flaton"), and the 

flaton supermultiplet X which is a singlet under the standard model gauge groups has the 
following superpotential: 



\ k Jf(n+3)fc 

W = ^{n + ?>)kM^ k -* +Cl (13) 



where A& denotes the coupling constant and we take Ai = 1. Here we impose the dis- 
crete Z n+ 3-symmetry to guarantee the flatness of the potential. In order to cancel out the 
cosmological constant the constant term C should satisfy 

\C\ ~ m 3/2 M*. (14) 



The cutoff scale of the model is denoted by M* in Eq. (pT3|) f| Note that we only have to 
consider the leading term of k = 1 since the higher power terms are highly suppressed. 
Then we obtain the following effective potential at low energy (\X\ <C Mq, M#) as 



n 



Va &) = v t -*tm> + — M — - — + 4^-, He) 

where we use the same letter X for the scalar component of the superfield. In this potential 
we have assumed the negative mass squared -Wq at the origin which is induced by the SUSY 
breaking effects and mo is of the order of the electroweak scale Aew ~ 100 GeV. Then one 
finds the vacuum expectation value (vev) of the flaton ((X) = M) as 



M 



n+l 



M? 2(n + 2) 



nm 3/2 + V n2m 3/2 + 4 ( n + 2 ) m o • (I 7 ) 



Therefore, the flaton is found to have a very large vev M ^> Aew since M* 3> Aew- Here 
notice that we assumed the soft SUSY breaking mass (mo) for the flaton, and hence the 
potential ([16|) is applicable only for the scale below masses of messenger fields when we 
consider the light gravitino (modulus) mass region in the GMSB models. Therefore, the 
vev of the flaton M should be smaller than their masses. However, in the original thermal 



9 As the cutoff scale of the theory it is natural to choose the gravitational scale 
term of the superpotential ( |T3| ) should be considered as 

1 V«+3 yn+3 

int ~ n + 3 M? ~ Mg ' 

Here A is introduced as a coupling which should be A < 0(1), thus Af* > Mq is naturally expected. 
Although we will take it as a free parameter in the following analysis, we will also mention about 
this lower bound on M#. 



. Thus the interaction 
(15) 
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inflation model, this constraint is always satisfied when we estimate the minimum of the 
modulus abundance in the next section. 

The vacuum energy of the flaton Vq is written as 



Vo 



n + 1 
n + 2' 



M 



ml + 



nm 3/2 



2(n + 2)(n + 3) 



nm 3 / 2 + \Jn 2 m\i 2 + A(n + 2)mfj 



(18) 



The masses of the scalar particles associated with X can be estimated by using the variables 
X and a which is defined as 



X 



M + 



X 

V2 / 



exp 



HI 



V2M 



(19) 



Here \ is the flaton particle which causes the thermal inflation. On the other hand, we 
call the imaginary part of X, a, as an "i?-axion". Since the superpotential ( |i~3"D posses 
an approximate U(1)r symmetry which explicitly breaks down to Z n+3 symmetry by the 
constant term C, a is considered as a NG boson with an explicit breaking mass proportional 
to C. The masses of the flaton and the i?-axion are estimated as 



m 2 = 2{n+ l)m 2 + 



n[n 



+ l)m 3 /2 



ni 



x ~ v " ' ^"" u ' 2(n + 2) 
2 _ n(n + 3)m 3/2 



nm 3 /2 + \j n2m l/2 + + 2 ) m o 



2(n + 2) 



nm 3 / 2 + \Jn 2 m\i 2 + 4(n + 2)ml 



(20) 
(21) 



Notice that if one takes m as 



ml > 



n 2 (3n + ll)(5n+ 13) 
4(n + l) 2 (n + 2) 



3/2; 



21 2 



for n = 1, 



(22) 



the flaton decay into two i?-axions is kinematically allowed. Therefore, in the GMSB models 
where the gravitino mass is in the range m 3 /2 ~ 10 eV-1 GeV, the decay process \ — > 2a 
seems almost to be open since m 3 / 2 m ~ A EW . On the other hand, in the HSSB models, 
we expect m ~ m 3 / 2 ~ A EW and the flaton might not decay into .R-axions. Although the 
decay channel of the flaton is crucial to estimate the dilution factor of the string modulus 
energy, we assume that the decay process \ — > 2a be always open in the whole gravitino 
mass region and that m always satisfy Eq. ( p2j ) in the original model. Furthermore, if 
mo ^> m 3 /2, the vev of the flaton M is determined independently on m 3 /2 as 



M 



1 



2(n + l) 



.n + 2, 

and we also obtain simple expressions for V , m x and m a as 

1 



Vn 



n 



ml ~ 



m 2 a ~ 



m'M z 
n + 2 

2(n + l)ml, 
n(n + 3) 



\fn + 2 



m 3/2 m . 



(23) 

(24) 
(25) 
(26) 
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We are now at the point to see how the flaton with the potential ([16]) causes a mini- 
inflation (= thermal inflation) at the late time of the evolution of the universe. For the 
thermal inflation to work, the flaton does not sit at the true minimum (|ITD, but sits around 
the origin due to the finite temperature effects in the early universe. To realize it, the 
flaton has to interact rapidly with fields in the thermal bath of the universe. The Yukawa 
interaction 

W = gtXg, (27) 

with £ and £ in the thermal bath is sufficient. Here g% is the coupling. When the flaton sits 
at the true minimum, the fields £ and £ obtain heavy masses ~ g^M. However, for the 
case that the flaton is near the origin, £ and £ are almost massless and can be in the thermal 
bath if they couple strongly to the fields in the thermal bath. This condition is satisfied if, 
for example, £ and £ are 5 and 5* in SU(5) gauge group respecting the unification of the 
gauge coupling constants. Then the interaction fl27|) gives an additional mass to the flaton 
in the early universe as 

n m i I V|2n+4 

Vp rr(X) = V + (cT 2 - ml) \X\ 2 - — ^ fx ra+3 + X* n+3 ) + , (28) 

where c is the order one coupling and we take c = 1 in the following analysis.^ Thus for 
the cosmic temperature T > T c (~ mo) the flaton sits at the origin. 

Here it should be noted that the Yukawa interaction (^7]) with g% ~ 1 is also required 



from another reason. In fact, the interaction (|27| ) with a large g% can induce the negative 
mass squared at the origin in the potential (|16|) by the renormalization group effects. 

During the time that the flaton is trapped at the origin, its vacuum energy Vq becomes 
dominant and the universe expands exponentially, i.e., the thermal inflation takes place 
fl8|| . There are three possibilities for the form of the energy which dominates the universe 
before the thermal inflation occurs. If the radiation energy dominates the universe before 
thermal inflation, the thermal inflation starts when Vq becomes comparable to the energy 
of the radiation at T = Tsti ~ Vq^- Then the thermal inflation lasts for Tsti > T > T c . 
On the other hand, if the universe before the thermal inflation is dominated by the energy 
of the modulus oscillation, the temperature at the beginning of the thermal inflation is 
estimated as Tsti ~ (Vq 7 '(w^g)) 1 ' 6 with the modulus mass m^. The universe before 
the thermal inflation might be dominated by the oscillating energy of the inflaton of the 
primordial inflation. In this case, Tsti ~ Vq^ 8 Mq A T]J^ where r vj is the decay rate of 
the inflaton (pi [see Eq. ( |B11| ) in Appendix |B|. Anyway, the flaton causes the thermal 



inflation for Tsti > T > T c and it occurs just before the electroweak phase transition since 
T c = m ~ A EW . 

We turn to the thermal history after the thermal inflation ends and see how the modulus 
abundance is diluted. When the cosmic temperature becomes smaller than T c , the flaton 
rolls down to its true minimum and oscillates around it. And the flaton decay occurs when 
the Hubble parameter becomes comparable to the total width T x of the flaton. By the 



10 In the case £ and £ are 5 and 5* in SU(5), the constant c takes a value c = 5g|/3. 
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flaton decay and the successive decay of the i?-axion, the vacuum energy of the thermal 
inflation is transferred into the thermal bath and the universe is reheated. At this epoch 
the tremendous entropy is produced and the string moduli is diluted significantly. 

Since we have assumed in the original thermal inflation model that m satisfy Eq. (f22[), 
the flaton always decays into i?-axions. This decay rate is given by 

where we have neglected the mass of the i?-axion. In addition, the flaton decay into two 
photons occurs through the Yukawa interaction fl27|) via one-loop diagrams of £ and £. For 
example, £ and £ are 5 and 5* in SU(5), its rate is estimated as 

a em \ 2 ml 



rtx - *T) = £ {^) j*. (30) 

Here we have neglected the effects of the SUSY breaking. In a similar way, the flaton can 
decay into two gluons if kinematically allowed, i.e., m x > 1 GeV, and the decay rate is given 
by 

Furthermore, the flaton might also decay into the gaugino pair. However, we forbid them 
since such a decay overproduces the LSPs and hence is cosmologically dangerous. 

On the other hand, the flaton has a chance to couple directly to the SUSY standard 
model particles. Within the renormalizable interactions the suitable charges of the Higgs 
supermultiplets H and H under the Z„ +3 -symmetry allow the superpotential 

W = X^XHH, (32) 

where A^ is a dimensionless coupling. In this case the ordinary u fi term" is forbidden but 
is induced by the flaton vev as \i = A M M, and hence the coupling constant A M should be 
extremely small so that [i becomes the electro weak scale q Through the interaction Eq. 
(|32|) the flaton can decay into Higgs bosons or Higgsinos if allowed. Their decay rate can be 
written as 

T( x -> 2h) = C h 1, (34) 

VA ; 16vrM 2 ' v ; 



The electroweak scale of /x can be naturally understood by the nonrenormalizable interaction [31] 



as 



(33) 

other than Eq. (]3^). However, the following discussion is almost same in both cases. 
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where we have ignored the masses of final states. Here Ch is a constant parameter and 
Ch — {^^M/ m x) A = (/V m x) 4 f° r the decay into the Higgs bosons and Ch = (A M /m x ) 2 
for the decay into Higgsino pair. In order that /i does not exceeds the electroweak scale 
A* ~ ^ew ~ m xi ^ should be Ch ^ 1. Since Higgsinos produced by the flaton decay might 
also lead to the overclosure of the LSP, we only consider the flaton decay into Higgs bosons. 
Note that when we take A M 7^ 0, the charged Higgs and Higgsino also contribute to the flaton 
decay into two photons and its rate becomes 

r(x-2 7) = i^) 2 4. (35) 

VA [; 72% V 4tt J M 2 K ' 

The .R-axion produced by the flaton decay has the decay processes similar to the flaton. 
The .R-axion can decay into two photons with the rate 

^em \ Win 



for A M = and 



r(a->2 7 ) = — ^ (36) 
V ' ; 9tt V 4tt J M 2 V ; 



for A M 7^ 0. In addition, it also decays into two gluons for m a > 1 GeV with the rate 

r(a-2,) = -Lf^) 2 ^|. (38) 
v y; 4tt M 2 v ; 

Note that the i?-axion can not decay into Higgs boson pair even for A M 7^ 0. Although it 
might decay into Higgsino pair, this process is assumed to be forbidden in the same way as 
the flaton decay. Therefore, the R-axion is assumed to have only the radiative decay modes. 

Now we are ready to estimate how an entropy is produced by the decays of the flaton and 
the i?-axion. When the Hubble parameter becomes comparable to the flaton's total width 
r x , the flaton decays into both the SM particles and the i?-axions. Since the i?-axion only 
has the interaction with the thermal bath suppressed by 1/M, the flaton energy transferred 
into the .R-axion could not reheat the universe at this time. Then the only the energy 
transferred into the SM particles reheats the universe at T = T S m by the flaton decay. The 
ratio of the entropy densities just before to after the flaton decay is estimated as 

A 5M = 1 + (1 - e a ) \ — Y° , (39) 

where e a denotes the branching ratio of the flaton decay into two i?-axions. 

After that, when the Hubble parameter becomes comparable to the total width of the R- 
axion, the .R-axion decays into the SM particles occurs and reheats the universe at T = T R . 
At this time the entropy of the universe increases by the factor A a : 

a 1+1 v ° ( 2m A 1 am 

a a 3 (2n 2 /±5)g*T?T R { m x ) A SM ' [ ^ } 
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Here note that the energy transferred into the i?-axions is diluted at the rate R~ 4 for T > m a 
while diluted at R~ 3 for T < m a . 

Then, all of the vacuum energy of the flaton Vq is released into the thermal bath and the 
thermal inflation does increase the entropy of the universe by the factor A: 

A = Asm x A a , 

= 1 + f1- 1 4 Vo 4 V ( 2m a \ 

+ 1 Caj 3 (27r 2 /45)^T|T5M ° 3 {2^/Ah)g*T*T R \m x J ' 1 ] 

Here the branching ratio of x ~ > 2a is e a ~ 1 for the case that the flaton can not decay 
into Higgs bosons, because the radiative decay channels of the flaton are only induced by 
one-loop diagrams and their rates are significantly suppressed. On the other hand, if the 
flaton decay into Higgs bosons is allowed, e a ~ 1/(1 + Even in this case, however, 

A takes its maximum value when e a ~ 1. This is because the i?-axion decay rate is much 
smaller than the flaton decay rate and hence m x T R /(2m a TsM) "C 1. Therefore in order 
to obtain the maximum entropy production, the flaton decay into Higgs bosons should be 
suppressed. In the following, in order to make a conservative analysis we take e a = 1 and 
use 

A = * * M • (42) 



3 (2tt 2 /45)^T3T r V m 



x 



For this (maximum) entropy production the reheating temperature Tr is determined by the 
total decay width of the i?-axion which can be written as 



3 



Cj±, (43) 



where the parameter C a depends on the decay mode and is given by [see Eqs. (|36"D , 
and (|38|)1 

, <M^) 2 forA, = 
Ca=\ 2 , (44) 



49 
72vi 



for the case that the i?-axion decays only into photons (m a < 1 GeV), and since a heavier 
i?-axion dominately decays into two gluons 

for m a ^ 1 GeV. Then, the reheating temperature Tr is obtained as 

T R ~ 0.96 vTVWg ^ 0.96 d' 2 ™* ™ G ■ (46) 
From Eqs. (E5), (p5|), and (|26|) we can write the vev and the vacuum energy of the flaton as 
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3/4^3/4^1/2 



, n 3/4( n + 3)3/4 m y* m »l* 

M ~ 96 C 1/2 - 1 j - 3/2 G (47) 

a (n + 2) 3 / 8 T fl ' 1 j 

° - °- 92 ^ ( W +V — if - ■ (48) 

The entropy production factor fl4*2|), therefore, is given by 

A ~ 2.0 x 10- "'("V'y' c. (^f (40, 

(n + 2) 2 \m / 

Note that A is independent on mo since mo — T c . It can be seen that the lowest reheating 
temperature Tr ~ 10 MeV gives the maximum entropy production as 

for the case n — 1 and the i?-axion decays dominately decays into gluons (m a > 1 GeV). 
Therefore, the thermal inflation is found to produce a tremendous entropy at late time of 
the universe and can dilute all of the unwanted particles which are long-lived like the string 
moduli. 

To end this section, we briefly discuss about the initial condition of the thermal inflation. 
In order to realize the thermal inflation, the flaton should be trapped at the origin of the 
potential by the thermal effects. Thus we have required the Yukawa interaction (p7[). If the 
flaton sits around the origin just after the primordial inflation, £ and £ become massless and 
give the flaton a mass comparable to T. Therefore the key point of the initial condition of 
the thermal inflation is where the flaton sits just after the primordial inflation. 

One may think that the flaton sits at the true minimum ((X) = M) just after the 
primordial inflation. In this case, since £ and £ obtain very heavy masses, the maxi- 
mum temperature Tmax achieved after the inflation should be higher than their mass, i.e., 
Tmax ~ — M in order to thermalize them. The maximum temperature is estimated 
as [see Appendix 

T MAX ~ 0.702 g^TMAx)- 1 ^ q^Trj) 1 ^ T^V/ /8 , (51) 

where the reheating temperature of the primordial inflation is Tri and its vacuum energy is 
Vj. Requiring Tmax ~ m ^ we obtain 

V} 1/4 > 2.03 g*{T MAX ) 1 ' 2 g*{TRj)- l ! A M 2 T^, 
*7 ™, ( M \ V Tri 



- 7 - 64Xl ° 7GeV lT0^GevJ ll^vJ • ™ 

Thus the inflation model with very low Tri should have very large vacuum energy. 

However, one can avoid this difficulty by considering the effects of the supergravity. 
The flaton can sit at the origin due to the additional SUSY breaking effect from the large 
vacuum energy of the inflaton. In this case the flaton starts to roll down toward its true 
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minimum when H ~ mo ~ ^ew- Then the flaton can be trapped at the origin, if the 
cosmic temperature at H ~ m is higher than the negative curvature at the origin, i.e., 
T(H ~ m ) ^ m . This gives the lower bound on the reheating temperature Tri as 

m 3/2 

Tri > 7.10 jjTT, 

m \ 3 / 2 



^ 459x10 GeV {mn^) ■ < 53) 

This condition is always hold since Trj should be higher than about 10 MeV from the 
BBN observations. Therefore, the initial condition of the thermal inflation can be explained 
naturally by the supergravity effects. 



IV. MODULI PROBLEM WITH ORIGINAL THERMAL INFLATION MODEL 

In this section, we estimate the modulus abundance in the presence of the original ther- 
mal inflation explained in the previous section and examine whether it could solve the 
cosmo logical moduli problem or not. 

If the Hubble parameter during the thermal inflation, -f/riQ, is larger than the modulus 
mass rricj,, the modulus oscillation starts after the end of the thermal inflation. The ratio of 
the energy density of this oscillation to the entropy density s after the reheating process of 
the thermal inflation is estimated as 

> 1.25 x HT 3 GeV ( ^f] ■ (54) 

where Tr is the final reheating temperature of the thermal inflation. Here we have used 
the fact that Tr > 10 MeV from the BBN observations and m x > 2m a since the flaton 
always decays into to the R-axions [see Eq. (^)j. We have also assumed that the modulus 
oscillation starts at least before the reheating of the thermal inflation, i.e., we have assumed 
the modulus mass as 

ms > l.l^f- ~ 4.5 x 10~ 23 GeV ( Tr ) . (55) 
* Mq V 10 MeV/ V ; 

The abundance (|54"|) becomes larger than the result Eq.(|4]) for Tri = Tr. Therefore, the 
cosmological moduli problem could not solved in this case and we should consider the case 

H TI < 77ty. 

In the following we will discuss the present modulus abundance in the presence of the 
original thermal inflation model with Hti < ttl^. Most of the previous works assumed that 
the oscillation of the modulus start after the reheating process of the primordial inflation 
completes. Thus the thermal inflation should occur after the end of the reheating of the 




12 The Hubble parameter during the thermal inflation Hti is estimated as Hti = Vq /(y/3Mo) 
with the vacuum energy of the flaton Vq. 
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primordial inflation. However, various models of the primordial inflation do not meet this 
assumption and give completely different consequences. In the present article, we will con- 
sider more general cases. The key time scales of the discussion are the following three scales: 
(T)H ~ : the cosmic time when the modulus starts to oscillate. (J1)H — the 
cosmic time when the reheating process of the primordial inflation ends. Here denotes 
the decay rate of the inflaton (pi of the primordial inflation. (Hl)iJ ~ Hti'- the cosmic 
time at which the thermal inflation occurs. Since we are considering the case > Hti as 
mentioned above, the thermal history of the universe is classified into the following three 
cases in general. 

• Case I: T Vl > > Hti- 

• Case II: > r w > H TI . 

• Case III: > H TI > Y ipi . 

For these three cases we will estimate the modulus abundance with the original thermal 
inflation model.Q 



A. Case I: For the case Y VI > > Hti 

First of all, we consider the case T VI > > Hti where a modulus field <^ starts to 
oscillate after the reheating process of the primordial inflation completes. In this case, the 
reheating temperature of the primordial inflation, Trj, which can be written as Eq. (^|), 
should be 



1/2 

T R i >Ts^ 7.2 x 10 8 GeV ' "'" > 



where denotes the cosmic temperature when the modulus starts to oscillate [see Eq. (Q)]. 
Therefore, a relatively higher reheating temperature is required and in some models of the 



inflation, such as chaotic or hybrid inflation, we can easily obtain such high Trj |29 . 

At T = the ratio between the energy density of the modulus oscillation to the entropy 
density is given by Eq. (j3|). After that, the energy density of the modulus is carried by the 
coherent oscillation. Hereafter, we call this modulus as "big-bang modulus" . In the presence 
of the original thermal inflation the present abundance^ of this big-bang modulus is diluted 
by the entropy production factor A [see Eq. (|42"D l and becomes 



13 The modulus oscillation is considered to start after the primordial inflation, i.e., < Hi = 

14 Here we also assume one modulus field with ~ m 3 / 2 to make a conservative analysis. 

15 For the unstable modulus of mass > 100 MeV, it corresponds to the abundance just before 
the modulus decay. 
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Moreover, it should be noted that the modulus energy is also produced after the thermal 
inflation. During the thermal inflation the modulus dose not sit at the true minimum but 
is displaced from it by an amount 5(p ~ (l^/m|M|l)</>o, and this causes the secondary 
oscillation of the modulus 181. We call this modulus as the "thermal- inflation" modulus. 



The abundance of the thermal-inflation modulus is given by 



PA 

s Jti 



IT/2 
2 V 



00 



X 



A' 



(58) 



Here notice that the thermal-inflation modulus starts to oscillate before the entropy produc- 
tion of the thermal inflation takes place and can be diluted. From Eq.([4~2"D the both modulus 
abundances are expressed as 



3.8- 



m 



qMq /T C V 



BB 



Tl 



0.38 



V 

V T R ( m 



\m J 




m\Ml \2m a 




Therefore, the total abundance of the modulus is given by 



^ -MAX 
s Jo 



BB 



TI. 



(59) 
(60) 

(61) 



We turn to estimate the lower bound on this total abundance of the modulus and compare 
with various cosmological constraints. The original thermal inflation model is parameterized 
by two mass scales mo and M* besides the gravitino mass m 3 / 2 (— m^). Here we take mo and 
Tr as two free parameters since M* is the function of mo, Tr and m.3/2. These two parameters 
are constrained as Eq. (^) for m since the flaton is assumed to decay into i?-axions in the 
original thermal inflation model, and Tr > 10 MeV from the BBN observations. 

The lower bound on the total abundance of the modulus (|H]) can be estimated by using 
the fact 




0.84 



S / BB \ S / TI 

(n + i)V2( n + 2)V4 m 2 TR / TO3/2 \ - - ( /; 



n 



V2(n + 3)V2 m 5/4 M 3/4 ^ 



V2 <T f /2 ( <Po V 



mo / \ M ( 



G , 



(62) 



where the equality holds when y^j — (v) > i- e -; 

1.9 + TrM 1 ^ /m 3/2 r 3/4 /^V /4 



m ° ~ Cl' 2 n 3 / 4 {n + l)V2( n + 3) 3 / 4 TO V 8 \ ) \ / 
Therefore, we obtain 



(63) 
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3 /™ ,\- 2 /T\3/i \ 2 



> 



2^9 {n + 2Y T A R m 3/2 \ / 



o ~ C a n 2 (n + l)V2( n + 3 )2 m ^ M 1 ' 2 \ J \mj \M, 



G, 



(64) 



and the lowest reheating temperature Tr ~ 10 MeV gives the lower bound on the total 
abundance of the modulus. It should be noted that the dependence of the index n in the 
flaton superpotential (|I3| ) appears only in the pre-factor. For n = 1 we write down it in 
terms of the density parameter as 

«^-»- ! (iiv)ia j © ! (!)' w 

for the case that a i?-axion can decay into gluons (m^ ^ 10 keV), and 



8.7 ■ 10'' (^r 3/ 7^V ^ 1^1 forA, = 



1 keV/ y y Vmo/ \M, 



A3 



'G , 



2.8 x 10 9 ( -P- ) [ ^ I I — ) I I for 



(66) 



1 keV/ V J Vmo/ \M, 



for the case that a i?-axion can only decay into photons (m^ < 10 keV). 

However, the above estimation should be changed when the modulus mass becomes 

nu > 1 2 TeV (n + 1)^ + 2)"* (*kV ^ 

^ ' n 14 /9( n + 3) 2 / 3 (3n + ll) 4 / 9 (5n + 13) 4 / 9 ^ ™> 4 ) \m ) ' 

-14/9 



/ \ -14/9 rp x 2/3 

~ 200 GeV ^ ( ± ) forn = 1, (67) 

\ J \m J 

because m in Eq. (|63| ) lies outside of the region (|22|) if one takes T R ~ 10 MeV. In such 
modulus mass region the reheating temperature which gives the lower bound on (j^j should 
be higher than 10 MeV and is given by 

n 7/4 (n + 3 )3/4 (3n + n )l/2( 5ri + 13 )l/2 m 9/8 / \ 7/4 ^ -3/4 



Ti? " 7 ' 1Xl0 (n + l) 1 /2( n + 2)iV8 Mi /8 V m W Vm 



and then we obtain from Eq. (p4|) 

Here note that the .R-axion is heavy enough to decay into gluons in the modulus mass region 

We show in Fig. [| the lower bound on the total modulus abundance. It is found that 
only the modulus with mass m^(~ m.3/2) ^ 100 GeV is cosmologically allowed. Thus the 
moduli problem in the HSSB models can be solved by the original thermal inflation model. 
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Original Thermal Inflation Model (Case I) 
n=l,b=l 




[GeV] 

FIG. 4. The lower bounds on the modulus abundance in the presence of the original thermal 
inflation for the case I: T Vl > > Hti- We assume = T c = mo and 4>q = Mq and take 

n = 1 and 6 = 1. The solid (dashed) line denotes the lower bound when A M = (A M ^ 0). Upper 
bounds from various cosmological constraints are all shown by the dotted lines. 



However, the light modulus ^ 1 GeV predicted by the GMSB models is still faced with 
serious cosmological difficulties even if one assumes the original thermal inflation in this case 
I: T VI >m^> H TI f23|. 

We have considered the parameters of the original thermal inflation model, m and M*, 
(i.e., m and T R ), as free parameters. Here we discuss how the lower bound on the modulus 
abundance obtained here is changed, when the cutoff scale of the original thermal inflation 
model is bounded from below as M* > M cr ~ Mq [see the footnote in Sec. ifj. This lower 
bound on M* leads to the lower bound on m as 



m > (n + 2) 1/2 



1.2 



C a 2 n 3 (n + 3) J 



n + l 



4(n + l) 
rp 3n-l 
J- D 



m 



3/2 



3(n + l) 2Q-1) 
3n — 1 3n — 1 



G 



Mr 



M, 



G 



(70) 



Since this lower bound becomes more stringent as n becomes large, which results in the 
larger abundance of the modulus, we take n — 1 in the following: 



m > 



3.2 x 10 



-2 



m| /2 V Mq 



(71) 



Due to this lower bound on m we find that in the modulus mass region 



ms < 0.18 GeV ( Tr ) 
VlOMeVy 



24/23 



-18/23 



T 

- 1 c 

m 



-6/23 



M, 
Mr 



16/23 



(72) 
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Original Thermal Inflation Model (Case I) 
n=l,b=l 
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FIG. 5. Same figure as Fig. || except for M* > Mq. We also show the results in Fig. |] by the 
thin lines. 



the abundance of the thermal-inflation modulus becomes always larger than that of the big- 
bang modulus. Therefore the minimum of the total abundance of the modulus is given by 
using Eq. (ffl]) as 

/ „ x A n w 1 n-7 T^15 /™ \ ~ n / n.r k 8 / ± \ 2 

(73) 



Pj \ ( PjL \ > 4.9 x 10- 7 Tft fm 3/2 \ 11 / M, 



Wo Wtj ~ C a 7 m^McV^W \MqJ \M g ) ' 
and the lowest reheating temperature Tr = 10 MeV leads to 

for m a > 1 GeV. On the other hand, for the heavier mass region, the minimum of the modulus 
abundance is same as the previous one Eqs.(|65D and (|69"D , even if one takes M* > M cr ~ Mq. 

We show the result in Fig. j|. One finds that the minimum of the modulus abundance 
becomes much larger than the previous results in the mass region given by Eq. ([T2|). On the 
other hand, the allowed region of the modulus mass ~ 100 GeV still survives even if one 
takes M* ^ M cr ~ Mq. Thus, for the case I: > > H^i the original thermal inflation 
can naturally solve the cosmological difficulties of the string moduli particles if the gravitino 
mass is m 3 / 2 (— m^) > 100 GeV which is predicted in the HSSB models, while cannot solve 
the moduli problem in the lighter gravitino mass region of the GMSB models. 
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B. Case II: For the case > T VI > Hti 

Next we consider the case that > F VI > H TI . In this case, the big-bang modulus 
starts to oscillate before the reheating process of the primordial inflation completes. Thus 
the energy density of the modulus is diluted by the primordial inflation as Eq. (|4]). Then 
the original thermal inflation takes place after that, and the mass density of the big-bang 
modulus is further reduced by the thermal inflation as 



with the reheating temperature 

T > T RI > 0.35V lf \ (76) 

where is the temperature when the big-bang modulus starts to oscillate and is given by 
Eq. @. This abundance takes its minimum value when T VI = H TI as 

J ±) > (&) = 4.4 x 10- V 1/4 ( 4f] 2 x i. (77) 

sJbB - \s) B Bm ° \M G ) A V ' 

On the other hand, the abundance of the thermal-inflation modulus is the same as Eq. ( |58|) 
in the previous case I. 

Let us estimate the lower bound on the total abundance of the modulus (|6lD in this case. 
One can write ( — ) and ( — ) in terms of mn and T R as 

V » J BBm V s ) T I U n 



2.2 (n + 2) 25 / 16 m 7 /8 T 5 R /2 ( %\ 



S J BBm Cl 1 ' n*/*(n + 1)V4(„ + 3)13/8 m 13/8 M 3/4 \ Mq 



(78) 



~024r »(^ + l) 3/2 (^ + 3) mgm 3 /2 fjQ_V ^ 
»/n ' (n + 2)3/2 m^MaTji \ Mq J [U) 

Both abundances become larger as m becomes larger, and ( — ) =( — ) is achieved 
when 

^ 2.01 (n + 2)«>/°° m^M^Tr 5 / ^\ 24/25 ^ 

° 1 ° Je9 " C a 14/25 n 2 V 2 5(n+l)l 4 /25( n + 3)21/25 m 21/25 { m J " W 

Therefore the minimum of the total modulus abundance is given by ( — ) with the lowest 

V s / BBm 

values of T R and m (< ( , m ) eq ). The lower bound on m Eq. (|22|) gives 



S / V S / BBm 

> _L2_ (w + 2) 9 / 8 (3n + ll) 7 / 16 (5n + 13) 7 / 16 T 5 R /2 f 7^\ 3 /^\ 2 

~ a 3/4 n 3/4 (n+1 )9/8 (n + 3)13/8 ^4^3/4 [ m J [ Mq ) 



(81) 
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For n = 1 the lowest reheating temperature Tr = 10 MeV leads to 



W> * 7, x 10- (^) 



~ 3 / 4 /m 3/2 V 3/4 /^N 3 



. nty j \m J \M G ) 
for the case that the i?-axion can decay into gluons (m^ > 100 MeV), and 

\ 2 



(82) 



n.h 2 > 



5.2 x 10 

2.3 x 10 



-2 ( m 4> 



V100 keV 
V 100 keV 



'm 3 / 2 ' 



-3/4 /... \ "3/4 x 3 / , x 2 



m J \M Gi 
T c V / _vo 



for A, 



(83) 



for 



G , 



for the case that the i?-axion cannot decay into gluons (m^ < 100 MeV). 

However, for the modulus mass region < 10~ 5 GeV, when the total abundance takes 
its minimum value Eq. (p3|), m becomes so small that the vacuum energy of the original 
thermal inflation Vq is less than the energy of the radiation at T = Tr. To avoid this failure, 
m should be 



m > 



1.5 



(n + 2) 1 ' 2 



T, 



12/7 



™3/2 



-3/7 



C 2/7 n 3/7( n + 1 )2/7( n + 3)3/7 m 3/7 M 2/7 ^ ^ 

Therefore the minimum of the total abundance is given by 

3.0 (n + 2) 2 T 4 /T^ 3 



(84) 



> 

o ~ C a n 2 {n + l) 1 /2( n + 3 )2 rn* /2 M G \m 



M, 



For n = 1 we obtain 



tt^h 2 > < 
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for A M = 



for A M ^ 



(85) 



(86) 



G, 



On the other hand, for the modulus mass region > 10 GeV, in order that m satisfies 
both mo < {mo)eq and Eq. (0), the reheating temperature should be higher than 10 MeV 



as 



T fi >7.7x 10- 4 



n 23/14 (n + 3 )3/4 (3n + H)25/56 (5n + 13) 25/56 ^15/14 
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Therefore we obtain from Eq. (|8T| ) 

.100 GeV J 



Oa/i 2 > 4.0 x 10~ 7 



m J 



M, 



for n = 1. 



g , 



(88) 
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FIG. 6. Same figure as Fig. |I| for the case II: > T Vl > Hti- 



We show the obtained lower bound on the total abundance of the modulus in Fig. |6]. It 
is found that since the primordial inflation does dilute the energy of the big-bang modulus, 
the lower bound becomes weaker than the previous case I and the allowed modulus mass 
regions are > 1 GeV and ~ 10 keV-1 MeV. Therefore, the gravitino mass region 
predicted by the HSSB scenario can be cosmologically viable and it should be noted that 
small window for the gravitino mass range in the GMSB models does appear. This feature 
is crucially different from the previous results. 

However, when we take the cutoff scale of the original thermal inflation model as M* > 
M cr ~ Me-, this new window is disappeared as shown in Fig. [F[ The lower bound on 
the total abundance of the modulus becomes more stringent for < 1 GeV where the 
abundance of the thermal-inflation modulus is always larger than the big-bang modulus one 
and then the limit is the same as Eq. (|T4| ) in the previous case. Therefore, in order that the 
original thermal inflation dilutes sufficiently the light modulus of mass ~ 10 keV-1 MeV 
in GMSB models, the extremely low cut off scale as M* ~ 10 6 -10 10 GeV is required. On 
the other hand, the allowed region for •> 1 GeV still exists even for the case M* •> Mq. 



C. Case III: For the case > Hti > r w 

Finally we consider the case > Hti > where the reheating process of the pri- 
mordial inflation completes after the thermal inflation ends and its reheating temperature 
becomes extremely low. In this case the present abundance of the big-bang modulus is given 
by [see Appendix 
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FIG. 7. Same figure as Fig. ^ except for M* > Mq. We also show the results in Fig. || by the 
thin lines. 
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This ratio takes its minimum value when T v = Hti as 
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As well as the previous cases we can rewrite in terms of mo and T R as 
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S J BBm 
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C a n 2 (n + iy/ 2 (n + 3) 2 m 2 /2 M G \m J 



M, 



G, 



(89) 



(90) 



(91) 



On the other hand the abundance of the thermal-inflation modulus is the same as the 
previous two cases [see Eq. (|7|)]. Now the condition (^) Bfim = (jf) 

holds when 



m = (mo) 



3.0 



(n + 2) 



7/6 



2/3^4/3 



ml T, 



Ca /3 n(n + l) 2 / 3 (n + 3) m 3/2 \m J 



(92) 



Therefore the abundance of the big-bang modulus (|91"| ) with the lowest values of T R and 
mo(< (mo)eq) gives the lower bound on the total modulus abundance flBTp. From the lower 
bound on mo [Eq. (E3)] in order to allow the flaton decay into .R-axions we find 
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For n — 1 the lowest reheating temperature Tr = 10 MeV gives 



Os/i 2 > 0.96 x 10" 



10 ( m <t> 



1 GeV 



™3/2 



mj VM, 



G , 



(94) 



when the i?-axion can decay into gluons (m^ ^ 100 MeV), and 

-i , m „ A / / \ 2 
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(95) 



when the i?-axion cannot decay into two gluons (m^ ^ 100 MeV). 

In the same way as the previous case II, however, for the modulus mass < 10 keV 
we find from the lower bound on mo [Eq. (|84|)1 
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and taking Tr = 10 MeV we obtain for n = 1 
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(97) 



When the total abundance takes its minimum value, the reheating temperature becomes 
higher than 10 MeV for > 3 GeV as 



T R > 7.0 x 10" 4 



3.4 GeV 



n 3 ' 2 (n + 3) 3 /\3n + ll) 3 / 8 (5n + 13) 3 / 8 ™ 3 % / T c \ - 1 
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(98) 



and from Eq. (^) we obtain for n — 1 



Oa/i 2 > 3.7 x 10" 



"fy \ / ™3/2 
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1 TeV/ \ 

where the flaton dominately decays into gluons. 



mo J \M, 
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FIG. 8. Same figure as Fig. Q for the case III: > Hti > 



We show the lower bound on the total abundance of the modulus in Fig. |8]. It is found 
that the lower bound becomes smaller than those in the previous two cases, because the 
big-bang modulus is also diluted by the primordial inflation during the thermal inflation. In 
this case the allowed regions for the modulus mass are ~ 1 keV-10 MeV and > 300 
MeV, and the allowed mass region in the GMSB models extends wider than the case II. 

In these allowed modulus mass regions, in order to obtain the lowest modulus abundance 
the required reheating temperature of the "primordial inflation" is Trj ~ 10 MeV for ~ 
1 keV-10 keV, and 



RI 
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1 MeVy 
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5/4 

! 

5/4 



(100) 



for m<^ ~ 10 keV-10 MeV. On the other hand, for the heavier modulus of mass > 3 
GeV, Trj is given by 



T, 



HI 



5.5 x 10 6 GeV ( -^-) 
VI TeV/ 



3/4 (rrn/X' 2 
\m<t> ) 



c 

m 



(101) 



Therefore, extremely low reheating temperature Trj ~ 10 MeV-10 GeV is required to dilute 
sufficiently the light modulus of mass ~ 1 keV-10 MeV, even if one assume the original 
thermal inflation model. 
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FIG. 9. Same figure as Fig. || except for M* > Mq. We also show the results in Fig. 
thin lines. 



by the 



Furthermore, if one takes the cutoff scale as M* > Mq, the allowed region for the 
modulus mass ~ 1 keV-10 MeV vanishes as shown in Fig. |9], and only the modulus 
of mass > 1 GeV is allowed. As well as the previous cases I and II, the lower bound 
on the total modulus abundance when M* > Mq is given by Eq. ( [71D for < 1 GeV. 
Therefore, in order to dilute sufficiently the light modulus, the thermal inflation model with 
the extremely low cut-off scale M* ~ 10 5 -10 12 GeV is required for ~ 1 keV-10 MeV, as 
well as the low reheating temperature of the primordial inflation as T&r ~ 10 MeV-10 GeV. 



V. MODIFIED THERMAL INFLATION MODEL 



As showed in the previous section, the original thermal inflation can dilute the relic 
abundance of the string modulus significantly, and the modulus mass regions m^i^ m 3 / 2 ) ~ 
1 keV-10 MeV, and ^ 300 MeV survive the various cosmological constraints if we 
consider the model of the primordial inflation with an extremely low reheating temperature. 
However, if we take the gravitational scale as the natural cutoff scale of the thermal inflation 
model, i.e. M* > Mq, the former allowed region vanishes and only the modulus mass region 
77i0 > 1 GeV predicted by the HSSB models is allowed. To dilute sufficiently the light 
modulus particle predicted by the GMSB models, one needs the relatively low cutoff scale 
M* ~ 10 5 GeV-10 12 GeV for ~ 1 keV-10 MeV. This weak point of the original model 
comes from the fact that the flaton can always decay into i?-axions. The vacuum energy of 
the thermal inflation is mostly transferred into relativistic i?-axions by the flaton decay and 
the i?-axions lose their energy faster than non-relativistic particles as the universe expands, 
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which leads to much less entropy production (i.e., less dilution of the modulus density). 
Therefore, one might expect to dilute the modulus density more effectively by the thermal 
inflation model which the decay process x ~ > 2a is not kinematically allowed. 

Furthermore, the original model is faced with another serious difficulty Since the poten- 
tial of the flaton possesses the exact discreet Z n+3 symmetry to ensure the flatness of the 
potential, the degenerate minima leads to the domain wall problem. In order to avoid this 
problem one has to introduce a term which breaks the symmetry explicitly in the potential. 

One of economical modifications of the original thermal inflation model, which solves 
above two difficulties simultaneously, has been proposed by Ref. |23|| . A linear term which 
breaks the Z n+3 completely is added to the original superpotential of the flaton Eq. (O) as 



W 



X 



n+3 



+ C + aX, 



(n + 3) M? 

with the dimensionfull parameter a which is required to be 

ml /2 m x M 

\ry\ ^ 



(102) 



(103) 



in order to eliminate the domain walls [j3~0|1 . Then this superpotential gives the low energy 
potential to the flaton as 



V eS (X) = Vo -^(X + X*)- ml\X\* + ^(X™+ 2 + X*™ +2 ) 



Mr 



n 



C 



n + 3 MqM™ 



+ 



\X 



2n+4 



M 2n 



(104) 



The vev of the flaton is estimated as 



(X) = M 



and the vacuum energy Vq is 



1 



2(n+l) 



(n + 2)(l-ar) 



n(l — x) + 1 
(n + 2)(l -x) 



(105) 



m 2 M 2 . 



(106) 



with the dimensionless parameter x < 1 which is defined as a = —xM n+2 /M™. Here we 
assumed m 3> m 3 / 2 . It should be noted that even in the presence of the explicit breaking 
term in Eq. ( |102| ) the dynamics of the thermal inflation is shown not to change much if 



m ^> m 3 / 2 and x is not so close to one |23 . 

In the flaton potential ( |lQ4j) we assume the SUSY breaking mass m at the origin as 
well as in the original model. Therefore, when we work in the light gravitino mass region of 
the GMSB models, we obtain the constraint that the vev of the flaton M should be smaller 
than the masses of the messenger multiplets. In the modified model, this constraint plays a 
significant role when we take the cutoff scale as M* > Mq [see the discussion in Sec. |yT[| . 

Then, masses of the flaton and the i?-axion are given by 
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Therefore, if x takes a value 



2{n + 1) _ 

(1 ^jX ^ — ; — — = X m i n j 
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(107) 
(108) 

(109) 



the flaton decay into two i?-axions is kinematically forbidden.PI Therefore, this modified 



inflation model with mo ^> 7773/2 and x in the region (109) gives one way to forbid the flaton 
decay into i?-axions and as well as to solve the domain wall problem. In the following we 
take x = x min for simplicity and m as 



„ n 2 (3n + ll)(5ra+ 13) a 
™o > \, , -^7 , ^ rni /2 



4(n + l) 2 (n + 2) 



t 3/2 



21 2 r 

— m 3/2 forn 



1. 



(110) 



for the purpose of the comparison with the previous results in the original model. This is 
the same condition for mo as assumed in the original thermal inflation model [Eq. (^2])]. 

Then we turn to discuss the thermal history after the modified thermal inflation ends. 
In the modified model the vacuum energy of the flaton is completely transferred into the 
thermal bath when the Hubble parameter becomes comparable to the total width of the 
flaton T x and reheats the universe at the temperature T = Tr. As well as the case of the 
i?-axion, the total width T x can be written as [see Sec. [I I ! 
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X M 2 ' 



111) 



Here C x is given by 



r 2 
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97T V 4:71 



for A, 
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' 1 forA^O 
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l 72tt V 4vr 



for m x < 1 GeV, 



1 /«,^ 2 



-V 



47r \AttJ 

for 1 GeV < m x < 2m h {m h is the Higgs boson mass.Q), and 



;ii3) 



16 In this region of x, the domain wall problem is solved, i.e., a satisfies Eq. ( |1Q3| ) . 
17 It should be regarded as the lightest Higgs boson mass and we take ~ 70 GeV. 
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for m x > 2mh. Then the reheating temperature is obtained as 



T R = 0.9dy/r x M G . (115) 

By using this reheating temperature, the modified thermal inflation model increases the 
entropy of the universe by a factor 

4 V n 

A~1H — (116) 



which is obtained by putting e a = and substituting T R for T RyS M in Eq. (J4TD . Comparing 
with the entropy production factor Eq. ( ^2[ ) in the original thermal inflation model you 
can see that the suppression factor {2m a /m x ) is dropped off so that the relic density of the 
string modulus is expected to be diluted more effectively. 



VI. MODULI PROBLEM WITH MODIFIED THERMAL INFLATION MODEL 

In this section we examine whether the modified thermal inflation model could solve the 
cosmological moduli problem or not. Crucial difference from the original model is that the 
flaton decay into two .R-axions is kinematically forbidden. Then the reheating temperature 
T R and also the entropy production factor A of the modified thermal inflation are determined 
by the flaton decay. 



A. Case I: For the case r v/ > ?ru > Hti 

First, we consider the case I: T lfi[ > > Hti- I n the presence of the modified thermal 
inflation model with the entropy production factor ( p.l6| ), the abundances of the big-bang 
modulus and the thermal inflation modulus are given from Eqs. (0) and ( j55[) as 



-3.8 " (p ' " "° R (—\ I ^j- 1 , (117) 
bb V Q \mj \ Mq J v 1 



V T R 



2 



-0.38^^1-^1 . (118) 
ti m\Ml \M G ) 

Then, as well as for the original thermal inflation model, we find the minimum of the total 
modulus abundance given by Eq. (|6lf) in the parameter space of m and T R . 

The vacuum energy of the modified thermal inflation ( |106j ) can be written in terms of 
m and T R as 



m 5 M G 

V0 — ^VQ 

Here CVo is defined as 



Cyo^P. (H9) 
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From Eq. ( |119| ) we can express both abundances in term of mo and Tr as 
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The dependence on tuq tells that the total abundance takes its minimum value when 
= (* 



is achieved, i.e. 
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(123) 



and we obtain 
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Therefore the lowest reheating temperature T R = 10 MeV gives for n = 1 
6.0 x 10~ 3 



n^h 2 > 



m^ 
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for A M ^0 



for m x < 1 GeV (the flaton decays only into two photons) and 

n^> 1.01x10- <J^-Y IU (^f l7 (*L 
* VI GeV J \mo) \M C 



(127) 



for 2mh > m x > 1 GeV (the flaton decays dominately into two gluons). 

On the other hand, when the flaton can decay into Higgs bosons (m x > 2m^\ one more 
free parameter Ch appears in the flaton decay width and the condition i^j) BB = (v) > 
i.e., Eq. (|123|) , could be achieved by taking a moderate Ch- Then we find from Eq. (|124|) 
that the lowest values of m which corresponds to m x = 2m^ and Tr = 10 MeV give the 
minimum of the total abundance as 



M, 



G . 
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FIG. 10. The lower bounds on the modulus abundance in the presence of the modified thermal 
inflation for the case I: T ipi > > Hti- We assume = m.3/2, T c = mo and </>o = Mq and take 
n = 1 and b = 1. The solid (dashed) line denotes the lower bound when A M = (A M ^ 0). Upper 
bounds from various cosmological constraints are all shown by the dotted lines. 

for n = 1. In fact, this result gives the absolute minimum abundance for the modulus with 
a mass ~ 6 - 26 GeV. (See Fig. [To] .) 

However, for the modulus mass region > 60 GeV, the reheating temperature Tr = 10 
MeV could not satisfy (ffj = 
Eq. (|m|) T R should be 

-3/4 /_ \ 7/4 

T R > 18 MeV 



because of the lower bound on m [Eq. ( |110| )1- From 



9/8 

\m J 



m 3/2 



100 GeV/ \m / \ 
Then from Eq. ( |124|) this gives the minimum abundance as 

\ 15 /8/T,\ 3/ Vm 3 / 2 
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100 GeV 
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for n = 1. 



g , 



(129) 



(130) 



Fig. [10] shows the lower bound on the total abundance of the modulus. We find that 
the modified thermal inflation can dilute extensively the modulus density than the original 
one in the light modulus mass region < 100 GeV. On the other hand, the lower bound 
is almost the same as the original model for > 100 GeVQ [see Fig. |[. We obtain 



18 The difference of the order one factor comes from the difference of the pre-factor in the expression 
of m x [see Eqs. (|^) and (|107|) 1 and so on. 
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FIG. 11. Same figure as Fig. 1C except for M* > Mq. We also show the results in Fig. [10| by 
the thin lines. The lower bound on the modulus mass [Eq. (|139[) ] is represented by the dot-dashed 
line. 



the allowed regions for the modulus mass: < 3 MeV and > 4 GeV. Therefore, the 
window for the gravitino mass predicted by the GMSB models does exist even for the case 
h r<^/ > m <t> > Hti and this point is crucially different from the original thermal inflation 
model. 

Let us discuss the effect of the lower bound on the cutoff scale M* > M cr ~ M G . In the 
same way as the original thermal inflation model, it leads to the lower bound on m : 

0.25 / M cr \ 1/2 

for n = 1. Here we consider only the case n = 1 since larger n induces more abundant mass 
density of the modulus. This lower bound on m does change the above analysis for the 
modulus mass region 

/ T \ 6 / 5 / T \ ~ 6 /5 , Tiir s 7/5 

m, < 4.9 x 10' 5 GeV ( ) ( ± ) ( ^ ) , (132) 

VIOMeV/ \m Q J \M G J ' v ; 



because the condition ( |123] ) could not be satisfied and be always larger than 

In this case we find from Eqs.( |122|) and ( |131|) 
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(133) 
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and the lowest reheating temperature Tr = 



10 MeV gives 



VI kevy vm g ; v m g/ 

for m x > 1 GeV. We show the lower bound on the modulus abundance in Fig. |TT] when 
M* > M G . 

Here we neglect the modulus mass region < 3 keV, since the vev of the flaton 
becomes always larger than masses of messenger fields to obtain the reheating temperature 
Tr ^ 10 MeV. In the GMSB models masses of messenger fields m mess can be written as 
m mess ~ (F mess )/A with the T-component vev (F mess ) in the messenger sector and A ~ 10 4 - 
10 5 GeV ||, and we obtain the upper bound on masses of messenger fields as 

m mess ~ ^ • (^35) 

On the other hand, the vev of the flaton can be written from Eqs. (|107|) and ( |115| ) as 

3/2 , ,1/2 

M = 3.40^ ™° _ G , (136) 
1 J- n 



for n = 1. Therefore, when one consider the GMSB models, in order that the vev of the 
flaton should be smaller than m mess we obtain 

0.64 m 2 ^M l J"TT n _ 
m ° ~ 7a7^ ~Wz • (137) 

X 

For the case M* > M cr ~ M G the upper bound on the reheating temperature can be obtained 
from Eq. ( |131| ) as 



Then, the reheating temperature Tr •> 10 MeV is achieved only when 



(138) 



A \ /M \ 3//4 

^ - -3/2 > 3.0 keV ( — — - , (139) 



10 4 GeV/ VM G 
where we consider the flaton dominately decays into gluons. 



From Fig. [TT| we see that the dilution of the modulus density becomes less effective for 



the modulus mass region Eq. (|132|) and the modulus with a mass < 3 keV is excluded. 



However, the allowed region of the modulus (gravitino) mass does exist in the region where 
the GMSB models predicts, even if we take M* > Mq in the modified thermal inflation 
model, while it does not exist in the original model even for the case III (m^ > Hti > T^) 
[see Figs. 0,0 and 0| . 
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B. Case II: For the case > T VI > Hti 

Next we consider the case > T VI > Hti- From Eq. ( [75] ) the modified thermal 
inflation model with the entropy production factor (|116| ) gives the abundance of the big- 
bang modulus as 

»\ =3 ,!»£ffi ! fAV (140) 
s)bb Vo \mj \M G ) 1 J 

On the other hand, the abundance of the thermal-inflation modulus is the same as Eqs. 
( |118|) and (|122|) in the previous case. Now the abundance of the big-bang modulus takes its 
minimum value when T VI = H^i as 

— ) > m .2.9^r^) 3 fAV (141) 

s)bb ~ \s)BBm Vq /A Vm / \M G J v ; 



Then the lower bound on the total abundance Eq. (|61|) is given by 



> 



s Jo - \ \sJ BBm \sJTi C'yf 3 mf 3 M G 8/23 W \M G 



with 



™» = ^r»^rQ) 12/23 . (Ms) 

Therefore, we obtain for n = 1 with the lowest reheating temperature Tr = 10 MeV: 

2 



O^ 2 > 



7.6 x 10-1 fZkf7# ! for A„ = 



10- 8 GeV 7 Vm o y V M g/ M 

4 . 3x io- 5( ^-)- 6/23 f^ 6o/23 fAy forW0 

V10- 8 GeV/ Vm o y V m g J 



(144) 



for m x < 1 GeV and 

for 2m ft > m x > 1 GeV. 

Furthermore, if the modulus mass becomes larger than about 10 GeV, the total abun- 
dance takes its minimum value when the reheating temperature is not 10 MeV but 



(146) 



due to the lower bound on m [Eq. flllOD I, which leads to 
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FIG. 12. Same figure as Fig. 1C for the case II: > F^j > Hti- 



n^h 2 > 1.6 x 10~ 6 



m4> \ 27 / 14 
100 GeV/ 



— ) 

m J 



6/7 



fn 3/2 



47/14 



M, 



for n = 1. 



G . 



(147) 



This result is the same as Eq. ( j88[) in the original thermal inflation model if we neglect the 
order one factor. We find that the dilution becomes ineffective when the flaton decays into 
Higgs bosons, and then the lower bounds listed above are the absolute ones. 

In Fig. O the lower bound on the total abundance is displayed. Comparing the result 
in the original thermal inflation model showed in Fig. |6|, the lower bound on the modulus 
abundance is extensively reduced for the lighter modulus mass region ^5 10 GeV, while it 
is almost the same for ^100 GeV. In particular, the lower bound becomes much weaker 
than the previous case I for the lighter modulus mass region < 1 GeV, so that only the 
modulus mass region m^(~ m 3 / 2 ) ~ 10 MeV-1 GeV is excluded by the stringent constraint 
from the cosmic X(7)-ray backgrounds. 

Moreover, when we take the cutoff scale as M* > M cr ~ Mq, the lower bound behaves 
as shown in Fig. This changes the lower bound for < 1 MeV and the bound is given 
by Eq. ( |134[ ) in the previous case I. Therefore, the modulus whose mass is m$ ^ 3 keV is 
excluded as well as the previous case I. 



C. Case III: For the case > Hti > 

Finally we consider the case > H TI > T VI , where the abundance of the big-bang 
modulus in the presence of the modified thermal inflation model is given by [see Appendix 
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FIG. 13. Same figure as Fig. 12 except for M* > Mq. We also show the results in Fig. 12 by 
the thin lines. The lower bound on the modulus mass [Eq. (|139[) ] is represented by the dot-dashed 
line. 



4.8 



BB 



T VI V£ ,2 M G WiJ 



M, 



This ratio takes its minimum value when V 



VI 



> 



BB 



S J BBm 



8.2 



= H TI as 
moT R f T C V 



V \m J \M Gi 



(148) 



(149) 



While the abundance of the thermal-inflation modulus is the same as the previous two cases 



(150) 



and given by Eqs. (|118|) and ( |122| ). Since i^j) BB becomes equal to (jy) TI when 



m 



C 



1-7 i /0 r 

~T/3 m <t> 2 R 
vo 



/3 T 2/3 ( Tc_\ 

m J 



2/3 



the lower bound on the total abundance ( pT| ) is given as 



> 



4.9 T, 



7/3 



C 2 JS mfM G 



rp . 10/3 



C 

m J 



Mi 



G, 



(151) 



Therefore, the lowest reheating temperature Tr = 10 MeV leads to for n = 1 

-1/3 / T \ 10 /3 / ' \ 2 



n.h 2 > 



3.4 x IO" 8 
1.6 x IO" 8 



1 keV 
1 keV 



-1/3 



T \ 
— ) 

rp . 10/3 

77lo/ 



M, 



G, 



G . 



for A M = 



for A M ^0 



(152) 
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FIG. 14. Same figure as Fig. 1C for the case III: ms > Hti > T 



for the flaton can only decay into photons (m^ < 1 MeV), and 
tt 6 h 2 > 1.6 x IO -11 



m„ 



1/3 ^,10/3 

100 Mo W \m ) \M G 



(153) 



for the flaton dominately decays into two gluons (m^ > 1 MeV). In the case III the flaton 
decay into Higgs bosons makes dilution ineffective and gives no absolute minimum of the 
total modulus abundance. 



We show the lower bound on the total modulus abundance in Fig. ITJ. It is seen that 
in the case III the modulus mass density is sufficiently diluted so that the whole modulus 
(gravitino) mass region 10 eV-10 TeV predicted by both models of the GMSB and the HSSB 
survives the various cosmological constraints. This is a crucial result when one considers the 
primordial inflation model with quite low reheating temperature. When the total abundance 
takes its minimum value, the required reheating temperature of the primordial inflation is 
estimated as, for the modulus mass region < 1 GeV, 



T, 



when the flaton decays only into two photons, and 

\ 5 / 12 /T r \ 5 " 



A M = 



(154) 



T, 



HI 



1.1 x 10 4 GeV 



100 MeV 



7 V 

— ) 

Too/ 



(155) 
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FIG. 15. Same figure as Fig. 14 except for M* > Mq. We also show the results in Fig. 14 by 
the thin lines. The lower bound on the modulus mass [Eq. (|139[) ] is represented by the dot-dashed 
line. 



when the flaton dominately decays into two gluons. On the other hand, for the mass region 
> 1 GeV, the required reheating temperature is given by 



m 



5.5 x 10 6 GeV 



\ 3 / 4 
1 TeV/ 



T 

± c 

m 



1/2 



fn 3/2 



1/2 



(156) 



where the flaton dominately decays into two gluons. Therefore, in order to sufficiently dilute 
the light modulus predicted by the GMSB models, the required reheating temperature of the 
primordial inflation is not so extremely low as that required in the original thermal inflation 
model [see Eq. fllOOD l. 



Furthermore, we show in Fig. [15] the lower bound on the modulus abundance for the 
case M* > M cr ~ Mq- The lower bound becomes more stringent for the modulus mass 

< 100 MeV where the bound is given by Eq. ( p8| ) which is the same as the previous two 
cases. However, we find that in the wide region of the modulus mass ~ 3 keV-10 TeV 
the cosmological moduli problem can be solved naturally by the modified thermal inflation, 
if the reheating temperature of the primordial inflation is low enough. 



VII. DISCUSSION 

The thermal inflation, whether it is original one or modified one, can dilute significantly 
unwanted long-lived particles (i.e., the string moduli). Especially, as shown in Subsec. |VT C| , 
the modified model can gives a solution to the cosmological moduli problem and the whole 
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modulus mass region m^(~ m.3/2) ~ 10 eV-10 4 GeV predicted by both the GMSB and 
HSSB models survives from various cosmological constraints [see Fig. [L4j], if the reheating 
temperature of the primordial inflation is low enough for its reheating process to finish after 
the thermal inflation ends. Furthermore, in this case, we find that even if one takes the 
gravitational scale as the cutoff scale of the modified thermal inflation model (M* > Mg), 
the modulus mass region m^(~ m 3 / 2 ) ~ 3 keV-10 4 GeV is allowed [see Fig. |15| . Note that 
gravitino is also diluted sufficiently and we have no gravitino problem. 

However, since the abundances of all relic particles are also diluted by the thermal 
inflation, one might be faced with the following problems; how generate the present observed 
baryon asymmetry and what is the dark matter of our universe. Here we give some possible 
solutions to them. 

We first discuss the baryon asymmetry of the universe. The generation of the observed 
baryon asymmetry Yg = (t2_b/s) ~ 10 _10 -10 -11 is one of the challenging question in particle 
cosmology. Furthermore, if one assume the thermal inflation in the history of the universe, 
the situation becomes worse. Because the primordial baryon asymmetry is also diluted after 
the thermal inflation and the temperature at that epoch is T ~ 10 MeV in order to dilute 
the moduli abundance maximally, the GUT baryogenesis and the electroweak baryogenesis 
do not work in this case. 



However, as pointed out by Ref. [19], enough baryon number could be produced by the 



Affleck-Dine mechanism |[L5|| . The generated asymmetry is related with the present modulus 



abundance as 19 



Y B <^(e±) = 3.6xl0-»(^)(-^)" ^ , (157) 
\ s J bb VI GeV/ \ m / 

where fi^ ss = p^ BB j ' p cr . Therefore, a tremendous baryon number is produced primordially 
(i.e., H ~ WI3/2) and enough asymmetry could be left even after the entropy production by 
the thermal inflation. Then the lower bound on the modulus abundance can be obtained 
from the present value of Y B as 

Y B \ ( \ ( m 3/2 



In Fig. [16] we show this bound together with the result obtained in Subsec. |VI C| [Fig 



rj|. You can see that this lower bound ( |158| ) lies above various upper bounds on the modulus 
abundance for > 4 MeV. Therefore the Affleck-Dine mechanism can generate enough 



baryon asymmetry in the lighter modulus mass region with < 4 MeV [12] even in the 
presence of the thermal inflation, and, on the other hand, we require other mechanism of 
baryogenesis in the heavier modulus mass region. 



In Ref. |3T| a variant type of the Affleck-Dine mechanism was proposed, where an .D-flat 
direction, LH U , could produce the required baryon asymmetry after the thermal inflation 
ends. They assumed that Trig + m 2 L be negative (m Hu and m L are the soft SUSY breaking 
mass of a H u and a scalar lepton doublet L) in order that the LH U condensate rolls away from 
the origin. However, in this case the potential along a specific direction is unbounded from 



below ||32[]. Therefore, this interesting mechanism of baryogenesis is not phenomenologically 
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FIG. 16. Same figure as Fig. 14. We also show by the dot-dashed line the lower bound on the 
modulus abundance from the present baryon asymmetry Yb = 1CP 11 which is generated by the 
Affleck-Dine mechanism. 



viable.0 

Although we have assumed the lowest reheating temperature Tr ~ 10 MeV of the thermal 
inflation to dilute the modulus density most effectively, one might have a chance to generate 
appropriate baryon asymmetry by the electroweak baryogenesis mechanism if one takes a 
reheating temperature as high as ~ 100 GeV. In Fig. [T7], we show the lower bound on the 
modulus abundance in the modified thermal inflation with > Hti > T VI and Tr > 100 
GeV. We find that the moduli with > 100 GeV is cosmologically allowed even if T R > 100 
GeV for A M = 0. Furthermore, for the case A M ^ 0, the allowed region becomes wider as 
> 5 GeV by the flaton decay into Higgs bosons. Therefore, in the heavy modulus 
(gravitino) mass region predicted by the HSSB models, one had a possibility to obtain the 
enough baryon asymmetry even in the presence of the thermal inflation, if the electroweak 
baryogenesis would work. 

Therefore, in the presence of the thermal inflation, we have two possibilities, so far, 
to produce enough baryon asymmetry as well as to dilute the mass density of the moduli 
sufficiently: (i) < 4 MeV by the Affleck-Dine mechanism and (ii) •> 5 GeV by the 
electroweak baryogenesis mechanism. 

Finally, we would like to mention about the dark matter candidates under the tremendous 
entropy production by the thermal inflation. Note that a dark matter of our universe is 
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We thank T. Yanagida for informing us of this point. 
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FIG. 17. The minimum abundance of the modulus in the presence of the modified thermal 
inflation for the case III: > Hti > The thick solid (dashed) line denotes the minimum 

abundance for the case = (A„ 7^ 0) when Tr > 100 GeV and the thin solid (dashed) line 
denotes that when Tr > 10 MeV. Upper bounds from various cosmological constraints are all 
shown by the dotted lines. 



diluted away as well as the string moduli. One possibility is the stable modulus itself, if its 
mass is less than about 100 keV, since in this mass region we can obtain ~ 1 without 
conflicting with the x(7)-ray background constraints. In this shown in Ref. [031 



the moduli dark matter with ~ 100 keV will be tested by the future x-ray background 
experiments with high energy resolution. The axion is another candidate for the dark matter. 
If its decay constant is high enough as J'pq ~ 10 15 -10 16 GeV, its energy becomes comparable 
to the whole energy of the present universe for the case that the reheating temperature is 
Tr ~ 10 MeV P|]]. Recently, the author and Yanagida proposed another candidate for 



the dark matter in the presence of the thermal inflation ||35|| . The superheavy particle of 
mass 10 12 -10 14 GeV, which was primordially in the thermal equilibrium, could be the dark 
matter, if its lifetime was longer than the age of the universe. 



ACKNOWLEDGEMENTS 

We would like to thank T. Yanagida for various suggestions and stimulating discus- 
sions. This work was partially supported by the Japan Society for the Promotion of Sci- 
ence (TA) and "Priority Area: Supersymmetry and Unified Theory of Elementary Particles 
(|{707)"(MK). 



44 



APPENDIX A: EVOLUTION OF ENERGY DENSITIES OF INFLATON AND 
RADIATION AFTER THE PRIMORDIAL INFLATION 



Here we briefly estimate the evolution of energy densities of the inflaton cpi and the 
radiation after the primordial inflation. The model of the primordial inflation is basically 
fixed by two parameters: the vacuum energy of the inflation Vj and the decay width of the 
inflation r w . 

When the inflation ends at t = tEi, the energy densities of the inflation p ipi and the 
radiation pr are given by 

P VI (t E i)=Vi, (Al) 
p R (t E i) ^ 0. (A2) 
The Hubble parameter of the universe is given by 

H{t " )sH » = & (A3) 

Then the inflaton causes a coherent oscillation for Hei > H(t) > F Vl . During this period 
the energy of the oscillation dominates the density of universe and the Hubble parameter is 
given as 

Hit) ~ P^ 1 ' 2 . (A4) 

Vsm g 

The evolution of p Vl and pr is governed by the following equations: 

p Vl + ZHp Vl = -T Vl p Vl , (A5) 

p Vs + 4:Hp VR = r vi p Vl , (A6) 

where the dot represents the derivative with respect to the cosmic time. When t ^> the 
solutions of them with the initial conditions ([Alp and are given by 



P.M^vJ^f] 3 , (A7) 



V R(t) J 
PR (t) ~ z/*r M G vj 



2 A, ^ T/1/2 (R{t E i) V /2 



'R(t E i)\ 5/ * 
, R(t) J 



(A8) 



where R(t) is the scale factor of the universe at the time t. Thus the radiation energy density 
for R(t) ^> R(tEi) is diluted by a rate R~ 3 ^ 2 (not R~ 4 ) as the universe expands. Note that 
Pr takes its maximum value 

p R \ MAX ~ 0.24T Vl M G Vj 1/2 , (A9) 

at R(t) ~ 1.5R(tEi)- Then the maximum temperature achieved after the primordial inflation 
is estimated as 

T MAX ~ 0.93 g.{T M j U c)- 1 ' A iy i A M%*V}'*, (A10) 

~ 0.70 g^TMAxY^g^TR^T^V,", (All) 
where Trj denotes the reheating temperature of the primordial inflation and is given by 



T RI ~ l.lg^TRjY^JY^MG. (A12) 
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APPENDIX B: ESTIMATION OF THE ABUNDANCE OF THE BIG-BANG 
MODULUS FOR THE CASE III: M > H TI > T VI 

Here we derive the abundance of the big-bang modulus for the case III: > H TI > T lf>1 , 
where the reheating process of the primordial inflation completes after the thermal inflation 
ends and its reheating temperature becomes extremely low. The mass density of the modulus 
oscillation is diluted not only by the thermal inflation, but also by the primordial inflation 
even during the thermal inflation. Therefore, we expect the abundance of the modulus is 
extensively reduced than case I and II. 

The big-bang modulus starts to oscillate with the initial amplitude O ~ M G at t = t BB 
when H ~ m^. At this time, the energy densities of the modulus, the inflaton, and the 
radiation are estimated as 



p Vl (t BB ) ^ 3m \M 2 Gl 
6 

Pr^bb) - -m^T^jM^. 





(Bl) 
(B2) 
(B3) 



Here we used the fact that the energy of the inflaton's oscillation dominates the universe for 
H E i > H(t) > and 



H{t BB f ~ m} ~ H 2 EI 



R{t 



EI 



Vi (R(t EI ) 



R(t 



BB 



3M G \R(t BB ) 



(B4) 



The cosmic temperature at t = t BB is estimated as 



Tbb - ^(Tbb)-/ 4 ^ 4 ^^^. 



(B5) 



The thermal inflation starts at t = tsri when the vacuum energy V of the flaton x 
begins to dominate the energy of the universe. Since we are considering the case > H TI , 
tsri should be tsri > t BB . Then, for t BB < t < tsri the energy of the universe is almost 
dominated by the inflaton energy as follows: 



P(t) = P<Pi(t) + PS) + PR(t), 
- P Vl (t)+ p<p(t), 



3mlM G + -m^M 2 G 



~ 3rr? M 2 ( R ( tBE 
- ^M G I 



( R(t B B) \ 3 
{ R(t) ) ' 



(B6) 



Therefore we obtain each energy density at t — t STI as 
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P x (t S Ti) = V (B7) 

M^)-^o(^) 2 (B8) 

p Vl (tsri)^V (B9) 

p R (t S Ti) ^ ^V^ /2 M G . (BIO) 



We find that the vacuum energy of the flaton when the thermal inflation starts is comparable 
to the energy of the inflaton's coherent oscillation. The cosmic temperature at t — tsri is 

Tsti - (^) 1/4 ^(T5T.)- 1/4 ry/\/ 1/8 M^. (Bll) 
While the thermal inflation lasts, the Hubble parameter takes the constant value given 

by 

T/l/2 

H{t) = H -=jk-a <B12) 
and the evolution of the energy densities are governed by the following equations 

p tf>1 + 3H T ip Vl = —T^jp^j, (B13) 
P(j} + 3H TlP<t> = 0, (B14) 
p R + 4:H TI p R = T^p^j . (B15) 

Here we have neglected the effect of the modulus decay. The solutions are 



R{tsTi) 



3 



^^{-W) • < B16 > 



r I ^ ( R(tsTl) 



R(t) 



(B18) 



The thermal inflation ends when the the cosmic temperature becomes T = T c . At this 
time (t = tETi) the energy densities are given by 

P X {tETl) = V Q , (B19) 

nit ^Vf 0o V Mi^S. (BOO) 

M - v ^Mfk> (B2i) 

= ^^( T C ) T c 4 - (B22) 
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Therefore, the ratio between the energy density of the big-bang modulus and the entropy 
density is estimated as 



T C V 



1/2 



8 %^T Vl M G \M G 



(B23) 



Below T = T c the flaton rolls down to its true minimum and oscillates around it. By 
the flaton decay (and by the i?-axion decay) the reheating process of the thermal inflation 
completes and the universe is finally reheated to the temperature T = Tr. Then the huge 
entropy is produced by a factor A given by Eqs. (|42|) and (116). Therefore, the present 
abundance of the big-bang modulus is given by 



TrVr 



c"0 



1/2 



BB 



8V3T Vl M G \M ( 



x 



L G . 



A' 



(B24) 



In the original thermal inflation model, the entropy production factor (f42|) leads to 



4.8- 



m^TR ( m v \ ( T r x ! 



BB 



V n 1/2 M G 



Vi v 



2m a 



m J 



M, 



G , 



(B25) 



On the other hand, in the modified model, Eq. ( |116|) gives 



4.8- 



BB 



V n 1/2 M G 



VI 



m J 



Mr 



(B26) 
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